GRADIENT ESTIMATES FOR PARABOLIC AND 
ELLIPTIC SYSTEMS FROM LINEAR LAMINATES 



HONGJIE DONG 

Abstract. We establish several gradient estimates for second-order di- 
vergence type parabolic and elliptic systems. The coefficients and data 
are assumed to be Holder or Dini continuous in the time variable and all 
but one spatial variables. This type of systems arises from the problems 
of linearly elastic laminates and composite materials. For the proof, we 
use Campanato's approach in a novel way. Non-divergence type equa- 
tions under a similar condition are also discussed. 
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1. Introduction 

We consider second-order divergence type parabolic and elliptic systems 
with coefficients and data which are irregular in one of spatial directions. 
This type of systems arises from the problems of linearly elastic laminates 
and composite materials. We are interested in the local regularity of the 
gradient of weak solutions to these systems. 

Problems of this kind have been studied by many authors; see, for in- 
stance, [31 [H HU [14] . In [3], Chipot, Kinderlehrer, and Vergara-Caffarelli 
considered the weak variational formulation of the equilibrium problem of 
a linear laminates, i.e. a domain in M. d consisting of a finite number M of 
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linearly elastic, homogeneous, parallel laminae. They proved that any weak 
solution u of the uniformly elliptic system div(v4Vu) = / is actually locally 
Lipschitz, under the conditions that / is in H k (Sl), k > [d/2], and the coef- 
ficients matrix A are constants in each parallel laminae. They also showed 
that the W 1,oc norm of u is independent of the number M, so that in the 
limiting case A are allowed to be functions of one direction alone. In |13j . 
Li and Vogelius considered scalar elliptic equations for a single real function 
u: 

div(AVu) = f + divg 

in a domain which consists a finite number M of disjoint sub-domains with 
C ,a boundaries. This equation models deformations in composite media 
such as fiber-reinforced materials. The matrix A and data are assumed to 
be Holder continuous up to the boundary in each sub-domains, but may have 
jump discontinuities across the boundaries of the sub-domains. Under these 
assumptions, the authors derived global W 1,co and piecewise C 1 ' estimates 
of the solution u for 5 G (0, d ^ +1 ^ \- Their results were later extended to 
elliptic systems for vector-valued function u by Li and Nirenberg |14j under 
the same conditions, and the range of 5 was also relaxed to (0, 2 ( a °fi) ]- 
The bounds obtained in [13, 14J, however, may depend on the number M. 
The corresponding problem for parabolic systems is more complicated. In 
a forthcoming paper, Li and Li |15j further extend some results in p3] to 
parabolic systems under an additional assumption that the coefficients and 
data are at least twice differentiable in t. See also [7] for a related result on 
parabolic systems. 

We would like to mention two recent papers [5] and |19j on "partial 
Schauder" estimates. In [5], Dong and Kim considered both divergence 
and non-divergence form second-order scalar elliptic and parabolic equa- 
tions. They proved that if the coefficients and data are Holder continuous 
in some directions, then derivatives of solutions in these directions are Holder 
continuous in the same directions. By using a different method, Tian and 
Wang [19] proved similar results for non-divergence form elliptic equations 
with coefficients and data Dini continuous in some variables. Under certain 
conditions, their results also extend to second-order fully nonlinear equa- 
tions. An interesting question is how much regularity one can expect in the 
"bad" directions. In this paper we address this question when the "bad" 
direction is one-dimensional. We note that the proofs in [5] and |19| do not 
apply to systems since the maximum principle is used in both papers. In 
the case that the coefficients are regular in all directions, a similar problem 
was studied long time ago by Fife [8] . 

In this paper, we are concerned with parabolic and elliptic systems: 



Vu := -ut + D a {A al3 Dpu) + D a (B a u) + B a D a u + Cu = div 5 + /, (1.1) 
Cu := D a {A^D p u) + D a (B a u) + B a D a u + Cu = div g + /. 
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The coefficients of V and C are assumed to be bounded, and the operators 
are uniformly nondegenerate. The aim of our paper is to obtain optimal 
regularity of weak solutions when coefficients are assumed to be regular in 
the time variable and all but one spatial variables. To be more precise, let 
us denote a typical point in M. d+1 by z = (t,x), where x = (x 1 , ■ ■ ■ ,x d ) := 
(x',x d ) and z' := (t,x r ). We shall prove that if the coefficients and data are 
Dini continuous in z' , then any weak solution u to (jl.ip is locally Lipschitz 
in all spatial variables, C 1 / 2 in t, and D x m and U := A^Dpu + B d u — 
are continuous; see Theorem 12,11 below for more precise statement. We also 
prove that a Holder regularity assumption in z' on the coefficients and data 
gives a better regularity of u. In particular, D x m and U are Holder in all 
variables; see Theorem 12.21 

In the special case that the domain consists of M parallel laminate sub- 
domains as in [3] mentioned above, we show that if the coefficients and data 
are regular in each sub-domain and may have jump discontinuous across the 
boundaries, then Du is regular in each sub-domain up to the boundary; see 
Remark 12.51 i). Thus our results generalize the aforementioned results in [3] 
by allowing more general coefficients and also deriving optimal C 1,<5 estimate 
for 5 € (0, 1). Unlike |13t 114] . we do not impose any restriction on <5, and 
the bounds of various norms are independent of M. However, it should be 
pointed out that although we allow subdomains to have curved boundaries 
(see Remark 12.51 ii)). the geometry of the domain considered in [13} 114] is 
more general than in the current paper. 

Our arguments are quite different from those in [13|. [T4"l [T5] and [5j [TU] . 
Let us give a brief description as follows. The proofs below are based on 
Campanato's approach, which was used previously, for instance, in [5| Hfi]. 
The main step of Campanato's approach is to show the mean oscillations of 
Du in balls vanish in certain order as the radii of balls go to zero. However, 
we are not able to follow this approach in the usual way due to the lack of 
regularity of u in the x^-direction. To overcome this difficulty, we appeal to a 
recent result in |4] regarding the L p estimate for systems with partially VMO 
coefficients. A crucial step in the proof is to deduce from this result some 
interior Holder regularity of D x m and U := A d @ Dpu for parabolic systems 
with coefficients depending on x d alone. We then use some perturbation 
arguments on D x m and U together with a certain decomposition of u to get 
the desired estimates. 

By using a similar idea, we also obtain the corresponding results for scalar 
non-divergence form parabolic equations: 

Pu := -u t + a ap D ap u + b a D a u + cu = f. 

We prove that if the coefficients and / are Dini continuous in z' , then any 
solution u of the above equation is C 1 in t, C ' in x, and i±t and D xx /u 
are continuous; see Theorem 16.11 Under the stronger condition that the 
coefficients and / are Holder continuous in z', we obtain additionally that 
ut and D xx >u are Holder continuous in all variables; see Theorem l6.21 These 
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theorems generalize some results in [HJ EJ HH] for the Poisson equation. In 
the case that the domain consists of M parallel laminate sub-domains and 
the coefficients and / are regular in each sub-domain, we show that the 
second derivative of u in the x rf -direction is also Holder continuous up to the 
boundary in each sub-domain. 

As mentioned above, some of our estimates rely on recent work about 
Lp-regularity for elliptic and parabolic equations (systems) with leading co- 
efficients VMO in some of the independent variables; see Section [3l This 
series of work was initiated by Krylov in [12J. For further developments in 
this direction, we also refer the reader to |1Q(, [J[ |6] and references therein. 

The organization of this paper is as follows. In Section [21 we state our 
main theorems for divergence form systems and introduce some notation. 
We prove some auxiliary estimates in Section[3j The proofs of main theorems 
are given in Section [5] and Section Finally we treat non-divergence form 
scalar equations in Section [6j 

2. Notation and main results 
We are concerned with parabolic systems 

Vu := -u t + D a (A aP Dpu) + D a (B a u) + B a D a u + Cu = div g + /, (2.1) 

where g = (<7i,<?2, - "' i9d)- The coefficients A a @ , B a , B a , C are n xn 
matrices, which are bounded by a positive constant K, and the leading 
coefficient matrices A a @ are uniformly elliptic with ellipticity constant v. 

for any ^ = (C*) G K dxn . Here 

u = (u\... g a = (gl--- ,gT, f = (f\ ■ ■ • , /T 

are (column) vector-valued functions. Throughout the paper, the summa- 
tion convention over repeated indices is used. We also consider the following 
elliptic system 

Cu := D a (A al3 Dpu) + D a (B a u) + B a D a u + Cu = div g + /. (2.2) 

In this case A a P , B a , B a , C, g, and / are independent of t and satisfy the 
same conditions as in the parabolic case. 

2.1. Notation. By Du = (D^u) and D 2 u = (Dyu) we mean the gradient 
and the Hessian matrix of u. On many occasions we need to take these 
objects relative to only part of variables. We also use the following notation: 

DfU = Uf, D x iU = U x r, D xx iU = u xx >. 

Set 

B' r (x') = {y G M.^ 1 : \x' - y'\ < r}, B r (x) = {y G R d : \x - y\ < r}, 
Q' r (t, x) = (t- r 2 ,t) x B' r {x), Q r (t, x) = (t- r 2 ,t) x B r {x), 
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and 

B' r = B' r (0), B r = B r (0), Q' r = Q' r (0,0), Q r = Q r (0,0). 

By N(d,p, • • • ) we mean that N is a constant depending only on the pre- 
scribed quantities d,p, • • • . For a (matrix-valued) function f(t,x) in R d+1 , 
we set 

(f)v = |p| jT /(*, x) dxdt= f(t, x) dx dt, 

where 2? is an open subset in R d+1 and |X>| is the d+l-dimensional Lebesgue 
measure of V. 

2.2. Lebesgue spaces. For p G (1, oo), we denote 

Wl' 2 {V) = {u: u, u t , Du, D 2 u G L P (V)} . 
We also denote HI" 1 (2?) to be the space consisting of all functions u satisfying 

inf {II#I|l p (£>) + INIl p (z>) \ u = dwg + h} < oo. 

It is easy to see that M~ l (V) is a Banach space. Naturally, for any u G 
M~ l (V), we define the norm 

IMIh- 1 ^) = inf { \\9\\l„(p) + \\Hl p {v) \ u = divg + h} . 
We also define 

nl(V) = {u: u,Due L p (V),u t G M^iV)} . 
We use the abbreviations Wp' 2 = Wp' 2 (R d+1 ) and U\ = ^(M d+1 ), etc. 

2.3. Partially VMO and partially Dini spaces. For a function u in V, 
we define its modulus of continuity u3 UjZ i (in the mean) with respect to z' by 

= sup I \Q r \~ 2 / / a/, z d ) - u(s,y',x d )\ 2 dydsdxdt J 

zoeR d+1 V JQr(z )nv JQ r (z )nv J 

r<R 

We say u is partially VMO with respect to z' if oj UyZ /(R) — >■ as R — > 0. 

We call a continuous increasing function uj on M+ a Dini function if cj(0) = 
and for any t > 

ft 

d(s)/sds < oc . 

'o 

We say function m in P is partially Dini continuous with respect to z' if its 
modulus of continuity u u ,z> in z' is a Dini function. In this case, we write 
u G C^ ini (P). Clearly, any function in C^ ini (X>) is partially VMO with 
respect to z' . 

In a similar way, in the time independent case we define oj U)X i and the 
space Cf} ni 



I[u](t) := f u(s)/sds < 
Jo 
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We note that our definition of Dini continuity is slightly different from 
the usual definition, where the modulus of continuity is measured in the 
uniform sense. 

2.4. Holder spaces. For 5 E (0, 1], we denote the C S / 2 ' S semi-norm by 

\u(t,x) - u{s,y)\ 

[U\5/2,S;T> ■■= SUp 

(t,x),(s,y)£V \t ~ S\ d < 2 + \x-y\ d 

and the C" 5 / 2 '* 5 norm by 

\ u \s/2,5;V '■= [ U ]s/2,5;V + \u\0;V, 

where \u\o-x> = sup-p \u\. Next we define 

N(l+5)/2,l+5;Z> '■= [ Du ]s/2,S;V + ( U )l+8;T> ' 
\ U \(1+S)/2,1+S;V '■= N(l+5)/2,l+5;D + |«|o;X> + \Du\o-T>, 

where 

/ \ \u(t, x) — u(s, x)\ 
Wi+s-.v : = SU P —r. Ki + ^/ 2 • 

By C7( 1+5 )/ 2 ' 1+5 (P) we denote the set of all bounded measurable functions 
u for which the derivatives Du are continuous and bounded in D and 

N(l+5)/2,l+5;D < °°- 

We define a partial Holder semi-norm with respect to z' as 

|u(t,a;) - u(s,y)| 

« *',*/2,*;I> : = SU P 77 1,5/2 , i f7) 

x d =y d ,(t,x)^(s,y) 

and the corresponding norm as 

\ u \z',6/2,S;T> '■= [ U }z',8/2,8;V + \ u \o;V- 

S /2 5 

By C 2 7 ' (P) we denote the set of all bounded measurable functions u on T> 
for which [u] z i^/2,5;T> < °°- Similarly, we define 

M 2 ',(l+j)/2,l+j;D : = [Ac'^z'.VV;© + ( U )l+8;V > 

and the space C^! +5 ^ 2 ' 1+5 (P). In the time-independent case, we also define 
Ha;',*; I ' \x',5 an d the space C^, in a similar fashion. 

2.5. Main results. We state the main results of the paper concerning di- 
vergence form parabolic systems. Roughly speaking, the first theorem reads 
if the coefficients and data are Dini continuous in z', then any weak solution 
u to (12.ip is Lipschitz in all spatial variables and 1/2- Holder in t. 
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Theorem 2.1. Let A G C°, im , B G Cj? m \ / G Loo(Ql) and g G C^ ni (Qi). 



Assume that u is a weak solution to (I2.ip in Qi. T/ien we /icwe u G 



C 1/2,1 (Qi/2) «nd 

Mi/2,l;Qi/ 2 ^ ^UH^'K 1 ) + l°lo;Qi + l/|o ; Q a + hlli 2 (Q 1 )), (2-3) 



where N = N(d, n, v, K, z' i^B z')- Moreover, D x m and U are continuous 



In the next theorem, we show that a partial Holder regularity assumption 
on the coefficients and data gives a better regularity of u. In particular, the 
spatial derivatives of u in x' and U are Holder continuous in all variables. 



Theorem 2.2. Let 5 G (0,1), A G C S J 2 ' S , B G C S J 2 ' S , f G Loo(Ql) and 



< N(\g\ g , mSjQl + |/|o ;Ql + ||u|U 2(Ql) ), (2.4) 
where N = N(d, n, 5, u, K, [A] Z , )S / 2>S , [B] z ',5/2,s)- 



Remark 2.3. It follows from the definition of U that D^u is continuous in 
z' in Theorem 12.11 and is in C J ' in Theorem 12.21 In other words, D^u is 
only discontinuous in the x d direction, which is a quite intuitive phenomenon 
from a physical point of view. 

Remark 2.4. The conditions on /, g and B in Theorems 12.11 and 12.21 can 
be relaxed. From the proofs below, it is easily seen that we only need / to 
be in some weaker Morrey space (cf. (|4.3[) ) . For g a ,B a ,a = 1, ••• ,d — 1, 
we only require the regularity with respect to the ^"-direction, in which the 
derivative is taken. 

Remark 2.5. i) Suppose that Qi is divided into M laminate sub-domains 
Di,T>2,--- ,T>m by M — 1 parallel hyperplanes with the common normal 
direction (0, • • • ,0,1). Under the conditions of Theorem 12.21 we assume 
in addition that A d/3 , f3 = 1, ••• ,d, B d and g d are in C 5/2,<5 (A) for each 
i = 1, 2, • • • , M, but may have jump discontinuities across these hyperplanes. 
Then we infer that 



is also in C s/2 ' 5 (T>i n Qi- £ ) for any z = 1, • • • , M and e G (0, 1/2), with the 

C S/2,S 

norms independent of M. Similarly, in Theorem 12 .11 if we additionally 
assume that A d P , /3 = 1, • • • , d, B d and are piecewise continuous, then D^u 
is also continuous in T>i fl Q\- £ - 



in Qx/2, where 



U := A dp Dpu + B d u - g d . 





8 



H. DONG 



ii) Our results can be extended to systems which model composite ma- 
terials similar to those considered in [13|. HI] . Let T > 0, and be a 
bounded domain in M. d which consists M disjoint subdomains Qi, ■ ■ ■ , Qm 
with C 1 ' 5 boundaries. Suppose that any point belongs to the boundaries of 
at most two of the subdomains. We also assume that A a/3 , B a and g a are 
in C S / 2 ' S ((-T, 0] x H~) for each i = 1, 2, • ■ ■ , M. For e > 0, denote 

tt e = { x G 0, : dist(x,<9f2) > e}. 

Let u be a weak solution to (|2.1|) in (— T, 0) x 17. Then by the standard 
technique of locally flattening the boundaries, we can apply Theorem 12.21 to 
obtain that 

u e C (1+5)/2 ' 1+5 ((-T + £,o] x (n £ nTk)) 

for any i = 1, • • • ,M and e G (0, 1/2). Similarly, if fi x , - - - , U M have C7 1,Dini 
boundaries and ^4 a/3 , i? a , g a are Dini continuous in (— T, 0] x for each 
z = 1, 2, • • • , M, then from Theorem 12. II we get 

u G C^'^C-T + e, 0] x (S2 £ n 0^)) ; » = 1, • • • , M, 

and Du is continuous in (— T+e, 0] x (f2 £ nf2j). However, the bounds of u may 
also depends on the distances between different inclusions. Therefore, this 
method does not apply to the case when the boundaries of three subdomains 
touch at some point. We remark that the estimates obtained in |131 [2] 
(see also recently [15] and [7]) are independent of the distances between 
inclusions, under a restriction on the range of 5. Estimates of this type are 
interesting from a physical point of view because, in the elliptic model of 
extreme valued conductivity, the gradient blows up as the distance goes to 
zero. 

Next we state the results for elliptic systems (|2,2p . which follow imme- 
diately from Theorems 12.11 and 12.21 by viewing the solutions to the elliptic 
systems as a steady state solution to the corresponding parabolic systems. 

Corollary 2.6. Let A G C$ im , B G C^ ni , f G Loo(-Bi) and g G C^ ni {B x ). 
Assume that u is a weak solution to (|2.2p in B\. Then we have 

\Du\ 0]Bl/2 < N(I[u}g iZ >](l) + \g\o;B! + \f\0;B! + \W\\ L 2 {B 1 )) , 

where N = N(d,n,v,K,u)A,z' ^B,z')- Moreover, D x m and U are continuous 
in B l/2 . 

Corollary 2.7. Let 5 G (0,1), A G C 5 X ,, B G C 5 X ,, f G L^Bi) and g G 
Ct.,(B\). Assume that u is a weak solution to (|2.2p in B\. Then we have 

\Du\ . Bl/2 + [D x/U ]s-b 1/2 + [U]s-,b 1/2 < N(\g\ x > iS;Bl + |/| 0; Bi + IMIl 2 (Bi))> 
where N = N(d,n,S, v, K, [A] x > jS , [B] X > :S ). 

Estimates in the same spirit as in Theorems 1 2 . 1 1 and 12 . 2 1 for non-divergence 
form equations are stated in Section [H 
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3. Some auxiliary estimates 
We will use the following property of Dini functions. 

Lemma 3.1. Suppose that uj is a Dini function, and 

oo 

w(t) := ak (u{b k t) XhH <i + u(l)XbH>i) , (3-1) 

fc=0 

for some constants a G (0, 1) and b > 1. Then oj is also a Dini function. 

Proof. Let 6j(t) = u(t) for t < 1 and &(t) = cj(1) for t > 1. Then Cj is 
bounded, uniformly continuous Dini function. Moreover, 

oo 

u(t) = ^2a k u(b k t). 

k=0 

Then it is easy to see that u is increasing and continuous. It follows from 
(ECU) that 

oo 
k=0 

for some constants N and 7 > depending only on a, b and wi. By Fubini's 
theorem, we also get /J oj(t)/tdt < 00. The lemma is proved. □ 

We will also need the following energy inequality and a variant of the 
parabolic Poincare inequality. 

Lemma 3.2. Let B = B = C = and R > 0. Suppose v € H^Qr) ^ s a 
weak solution to the equation 

jVv = divg + f inQ R ; 
\ v = on dpQn, 

where f,g S -^(Qi?)- Then we have 

\\ Dv h 2 (Q R ) +^ -1 |MIl 2 (q b ) < n \\9\\l 2 (q r ) +NR\\f\\ L2{QR) , 

where N = N(d, n, v) > 0. 

Lemma 3.3. Let p G (I, 00), r G (0, 00) and u G C% c (R d+1 ). 

i) Suppose B = B = C = 0, g,f€ L p j oc , and Vu = div g + f in Q r . 
Then 

[ \u(t,x)-(u) Qr \ p dz<Nr p [ (\Du\ p + \g\ p + r p \f\ p )dz, (3.2) 
where N = N(d, v,p) > 0. 
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ii) There is a constant N = N(d,p) such that 

j \Du(t,x)-(Du) Qr \ p dz<Nr p f (\D 2 u\ p + \u t \ p ) dz , 

\u(t,x) - (u) Qr - X P{D p u) Qr \ p dz < Nr 2p [ {\D 2 u\ p + \u t \ p ) dz. 

JQr 

(3.3) 

Proof. See, for instance, Lemma 3.1 and Lemma 3.2 of [12J. We remark that 
by an approximation argument, (|3,2p remains valid for u £ Hz loc , and (|3.3p 

1 2 

remains valid for u G W p \ oc . □ 

In the remaining part of the section, we shall prove some local estimates, 
which are deduced from the results obtained in [3J. 

Lemma 3.4. Let p G (l,oo). Assume A a/3 are partially VMO in z' , u € 
Hp(Qi) and 

Vu = divg + f, (3.4) 

in Q\, where f,g£ L p (Q±). Then there exists a constant N = N(d, n, u, K, uja.z 
such that 

!MI^(Qi/2) ^ N (W u h P (Q 1 ) + llsll-MQi) + ll/IU„(Qi))- 

Proof. The lemma follows from Theorem 2.2 of [4] by a standard localization 
argument. For the sake of completeness, we give a proof in the Appendix. 

□ 

By using the Sobolev imbedding theorem and a bootstrap argument, we 
get 

Corollary 3.5. Let p,q E (l,oo). Assume A a/3 are partially VMO in z' , 
u € Cfa. satisfies (|3.4p in Q\, where f,g G L q (Qi). Then there exists a 
constant N = N(d, n, v, K, oja,z' > Pi q) such that 

IMki(Q 1/2 ) < N(\\u\\ Lp{Ql) + lblU g (Qi) + ll/IU,(Qi))- 
In particular, if q> d + 2, it holds that 

\ u U/2,r,Qi/2 ^ N (W u h P (Qi) + IMU,(Qi) + ll/IU,(Oi))> 
where 7 = 1 — (d + 2)/q. 

Next we consider systems with coefficients depending on x d alone. We 
denote 

V u = -u t + D a (A a P(x d )Dpu), 

and 

U := A^Dpu, i.e., IP = A^D p u j , i = 1, • • • ,n. 
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Lemma 3.6. Let p G (l,oo). Assume u G C?° satisfies Vqu = in Q\. 
Then for any nonnegative integers i,j such that i+j > 1 and any q G (1, oo), 
there exists a constant N = N(d,n,p,q,i,j,v) such that 

\\ D HM\n(Q 1/a )<mDu\\ Lpm . (3.5) 

For any 7 G (0, 1), we also have 

\DlDiMy/2, r ,Q 1/2 < m\Du\\ Lp(Ql) . (3.6) 

Moreover, 

\\DiDi>U\\Hl(Q 1/2 ) < N\\Du\\ Lp{Ql) , (3.7) 

\DiL&U\ l/2tTtQl/3 < N\\Du\\ Lp{Ql} . (3.8) 

Proof. Note that Vq{D\D x iu) = in Q\. Thanks to Lemma 13.41 and Corol- 
lary [331 to prove (|3.5p it suffices to show that for any 1/2 < r < R < 1, 

lklU 2 (Q r ) < iV(d,n,i/,r,i2)||Z)u|| £ , 2 (Q Jl) . (3.9) 

Indeed, if i > 1, by Lemma 13.41 and (j3.9|) . for any 1/2 < n < ri < r3 < 1, 

IIA^>ll^(o n) < iV||^i}>|U 2(Qr2) < ^ll^r 1 ^.^!!^^)- 

Similarly, if j > 1, we have 

\\V\V 3 M\u\ { Q T1 ) < N\\DlDi lU \\L 2 (Q r2 ) < N\\DlDi7 l u\\ nl{Qrs y 
Repeating this procedure to reduce i and j, we reach 

\\DlD J x ,u\\ n i,2 {Qi/2) < N\\u t \\ L2{ Q 3/4) + N\\Du\\ L2{Q3/4) < N\\Du\\ L2(Ql y 



In the last inequality, we used f|3.9[) . For the proof of (|3.9p . see for instance 
Lemma 3.3 of [1]. Inequality (|3.6p is deduced from (13. 5p by the parabolic 
Sobolev imbedding theorem. 

To prove (|3.7p and (j3.8l) . one only needs to observe that in Qi, 



d-l d d-l d 

a=l/3=l a=l/3=l 

and 

d d 

D X ,U = A d ^D x ,Dpu, D t U = Y A^DtDpu. 

/3=1 /3=1 

The L p norms of the right-hand sides are bounded by ||Z>u||£ p due to (|3.5[) . 

□ 
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4. Systems with partially Dini coefficients 

This section is devoted to the proof of Theorem 12. 11 The following lemma 
reduces the estimate of [u] 1/2,1 to the estimate of |Du|o- 

Lemma 4.1. Let u be a weak solution to (I2.ip in Q\. Suppose \u\q-q 1/2 < 00 
and \Du\o-q 1/2 < 00. Then we have 

Ml/2,l;Qi/ 4 ^ N (\ u \0;Q 1/2 + \ Du W,Q 1/2 + \f\o-,Q 1/2 + Mo ; Q 1/2 )- (4-1) 
Proof. Rewrite (|2.ip as 

-u t + D a (A a PDpu) = D a (g a - B a u) + f - B a D a u - Cu := D a3a + f, 
where 

Qa = g a - B a u, f = f- B a D a u - Cu. (4.2) 
We fix zq G Q1/4 and take r € (0, 1/4). By Lemma we have 

f \u-(u) Qr{zo) \ 2 dz<Nr 2 [ (|Dni 2 +[g| 2 +r 2 |f| 2 ) 

JQr{z ) JQr(zo) 

< Nr n+4 (\Du\ . Ql/2 + Mo;Q 1/2 + \g\0;Q 1/2 + l/|0;Q 1/2 ) • 

The inequality (|4,ip then follows from Campanato's characterization of 
Holder continuous functions; see, for instance, [17\ Lemma 4.3]. □ 

Now we are ready to prove Theorem 12.11 We shall divide the proof into 
two steps. 

Step 1: Estimate of \Du\q. We first assume that all the coefficients and 
data are smooth, so that by the classical theory \Du\o ] q 3/4 < 00. We take 
< 71 < 7 < 1. Fix a point z G Q3/4, and take < r < R < (3/4- |x |)/2. 
Now take z[ G Q'r(z' ) and denote 

V z ,u := -u t + D a (A a P(z[,x d )Dpu). 

Recall the definitions of g and f in (|4.2p . Then we have 

V z '_u = div(g + m) + f, 

where 

ml l (z) = (Af(z[,x d )-Atf(z))D^. 

Let 

u (x d ) = J (A dd (z / 1 ,s))- 1 Q d (z' ll s)ds, u e = u-u . 

Clearly, 

V z <u e = div(g(z) -g{z[,x d ) + m(z)) + f. 
Let v be a weak solution to the equation 

V Z 'V = div(g(z) - q(z' , x d ) + m(z)) + f in Q R (z ); 

v = on d p Q R (z ). 
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By Lemma 13.21 we get 
\\Dv\\l 2 (q 

< N\\g(z) - 5{z[,x d ) + m\\ L2(QR(zo)) + NR\\f \\l 2 (q r (z )) 

< N\\q(z) -Q(z[,x d ) + m\\ L2{QR{zQ)) + N\^, QR{zo) R d l 2+1+ ^. (4.3) 
Let w = u e — v. Then w satisfies V^w = in Qr(zq). Denote 

W := A d P(z[,x d )Dpw. (4.4) 
It follows from Lemma 13.61 and a suitable scaling that 

/ \D x ,w - (D x ,w) Qr{zo) \ 2 + \W- (W) Qr(zo) \ 2 dz 

JQr(z ) 

< N{r/R) d+2+2 ~< f \Dw\ 2 dz 

< N{r/R) d+2+2 ~< I \D x ,w\ 2 + \W\ 2 dz. (4.5) 



In the last inequality, we also used the nondegeneracy ofA dd . Define 




and 

d-l 

w:=w~Y^ xP { D P w )q r {z ) ~ Kx d )- 
(3=1 

We define W and U e as in (|4.4p with w and u e in place of w respectively. 
Then, 

D x ,w = D x ,w - (D x ,w) Qr{zq) , W = W - (W) Qr{zo) . 
A direct calculation shows that w also satisfies V z >w = in Qr(z ). We 
substitute w and W in (|4.5p by w and W to get 

f \D x ,w - (D x ,w) Qr(zo) \ 2 + \W- (W) Qr(zo) \ 2 dz 

< N(r/R) d+2+2 ^ [ \D x ,w - (D x ,w) Qr(zo) \ 2 + \W - (W) Qr[zo) \ 2 dz. 

(4.6) 

We combine (|4.3|) with (|4.6|) and use the triangle inequality to obtain 




\D x ,u e - {D x ,u e ) Qr{zo) \ 2 + \U e - (U e ) Qr(zo) \ 2 dz 



< N^r/Rf^ [ \D x ,u e - {D x ,u e ) QR{zo) \ 2 + \U e - {U e ) Q[l{za) \ 2 dz 
JQr(zo) 

+ N\\g(z) - Q (z[,x d ) + m\\l (QR{zo)) + JV|f \l, QR{zo) R d+2+2 ^ • (4.7) 
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where N\ = N\(d,n, u). It is not very convenient to use (j4.7j) since both u e 
and U e depend on zq. However, by the definition of uq, u e and U e , 

D x m e = D x ,u, U e = A d P{z[,x d )Dpu - g d (z[,x d ) 

Recall that 

U = A d ^Dpu- Sd (z), (4.8) 
which is independent of the choice of Zq. Clearly, in Qr{zq) 

\U(z) - U e {z)\ < N\A{z) - Aiz^x^WDulo.Q^ + \g d (z) - g d (z[, x d )\. 

Thus, coming back to (|4.7j) we get 

\D x ,u - (D x ,u) Qr{zo) \ 2 + \U- {U) Qr{zo) \ 2 dz 

< N l {r/R) d+2+ ^ [ \D x ,u - (D x ,u) Qr(zo) \ 2 + \U- (U) Qr(zo) \ 2 dz 
JQr(zo) 

+ N U(z) - s(4,x d )\\l 2{QH{zo)) + N\\A(z) - M^IIUqhWWIqhM 
+ N \f\l,Q R ( Z0 ) Rd+2+271 - (4-9) 

Now we take the average of both sides of (j4.9|) with respect to z[ £ Q' R (zo), 
and use the definitions of g and f to obtain 

\D x ,u - {D x ,u) Qr{zo) \ 2 + \U- (U) Qr(zo) \ 2 dz 

.■(z ) 

< N^r/R)^ 2 ^ [ \D x ,u - (D x ,u) Qr{zo) \ 2 + \U- (U) Qr(zo) \ 2 dz 
+ N(u g , z ,(R) +u B ,z>{R)\uW, Qz/A +cj A)Z ,(R)\Du\ . Qr(zo ^ R d+2 

+ N (N 7i/2j7i;Q3/4 + |f| ;Q B ( 20 )) 2 ^ +2+271 - (4-10) 

Set r = tR for some r € (0, 1) to be chosen later, and denote 



0rOo) = + \D x ,u- (D x ,u) Qr{zo) \ + \U - (U) Qr{zo) \ dz. 

It follows from (|4TT0|) that 

<f> TR (zo) < N^Mzo) + N (oj 2 >zl (R) + Mo ; Q 3/4 ) 

+Nu%AR)\Du\l QR{zo) T- d - 2 + N (N 7i/2j7i;Q3/4 + \f\l )QR{zo) ) r- d - 2 R 2 ^. 

We fix r = r{d, n, v) < 1 sufficiently small such that N\t 21 < 1/2. By an 
iteration, we obtain, 

k 

<P TkR (z ) < 2~ k Mz°) + NY,^(^lAr k ' ] R) + " 2 B ,AT k - j R)\<,Q 3/4 

3=1 
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k 

+u;i^(^-^)|^|2 ;Qi?( ^ ) ) + Ar([u]^ l/2m;Q3/4 + |f| 0;Q ^ ( , 0) ) 2 ^2^(^-^) 2 ^. 

3=1 

(4.11) 

We define 

Vv(^o) = f \D x 'u\ + \U\ dz. 

J Qr(zo) 

By the triangle inequality and Holder's inequality, 



\tprr(zo) -1p r (z )\ < + \D X ,U - (D x ,u) QrM \ + \U - (U) Qr{zo) \ dz 

< N(Mzo)) 1/2 - (4.12) 
Combining (|4.11|) and (|4.12p . we deduce 

\^rM-^rM\ < N2- k / 2 <p){ 2 (z ) 

k 

+NY,2~ j/2 (u g ,z>(T k - J R)+U JB , z ,(T k ^R)\u\ 0;Q3 ^ 



3=1 

k 

+ n (n 7i/2 , 7i; q 3/4 + \f\o; QRM ) E 2-^(T*-^r. 

3=1 

Summing the inequality above in gives 

Vv*i?(^o) < ^(z )+A r 4 /2 (^o)+A r E (u g , z/ (r k - j R) + u> BtZ >(T k -iR)\u\ ., Q3/4 

3=1 

fe 

+iv£ c^r^'i^lo^o) + 2V (M 7i/2 , 7i; q 3/4 + |f|o;Q«(* )) 

3=1 

where iV = N(d,n, v). To estimate the summations on the right-hand side, 
we recall 

I[w](r) = / oj(s)/sds. 
Jo 

It is easy to see that for any Dini function u and k > 1, 

fe 

^w(t*-^) < NIlui^T^R). 

3=1 

Therefore, we deduce 

^ T * fl (z ) < V'i?(^o)+A r 4 /2 (^)+A r /K,'](r- 1 J R)+iV/[a; B , z ,](T- 1 ^)| U |o ; Q3 /4 

+ iVJ[a;^](r- 1 iJ)|£)«| . Ojl(a0) + N (M 7i/2 , 7i; q 3/4 + lf]o ; Q«(^)) R 11 ■ 

(4.13) 
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It follows from (I4.13|) and the definition of f that 

(\ Du \)Q rkR (zo) < Mza) + N4>]( 2 (z ) +NI[u; g , z ,](T- 1 R) 

+ WKz'K^^jHo^^ + N(\g\ 0] Q 3/4 + |n| 0; Q 3/4 + [u] 7i/2j7i; q 3/4 ) 

+ N^A^ir^R) + R^)\Du\ , Qr{zo) + N\f\ . QR(zo) RK, (4.14) 

where N 2 = -/V 2 ((i, n, i/) > 0. We choose Rq = Ro(d,n,u,ujA,z') £ (0> r /4)) 
such that 

N 2 (I[coa,z'](t- 1 R ) + R^ 1 ) < 2- d ~ 2 . 

Now we confine zq to := Q 3 / 4: _2-e-if io , i = 1, 2, ... and set i? = 2~ e ~ 2 Ro. 

Since and zo € in (|4.14p are arbitrary, due to the Lebesgue lemma, 
we obtain 

\Du\ . qW < AT 3 2^+ 2 )/ 2 (\\Du\\ L2iQl) + | 5 | 0; q 3/4 + I[u g>z/ ](1) + |n|o ;Q3/4 ) 

+ iV 3 2^ +2)/2 (|/|o;Qh(.o) + M7iA7i;Q3/ 4 ) + 2- d - 2 \Du\ 0]Qie+lh (4.15) 

where N$ = N^(d, n,u, 0Ja,z'^B,z') > 0. Multiplying both sides of (|4.15p by 
2~e(d+2) anc j summ i n g j n l = 12, ... yield 

CO 

^2-^+ 2 )| J D U | 0;QW < iV 3 (\\Du\\ L2(Ql) + |g| . Q3/4 + l«|o;0 3/4 ) 

CO 

+ iV 3 (|/|o ; Qh(, ) + M-nAwQ./*) +E2~ (m)(d+3) l^«lo i 0(«-i)- 

i=i 

Since we assume |_Du|o ; q 3/4 < oo, it follows from the inequality above by 
absorbing the summation on the right-hand side into the left-hand side that 

\Du\ . QW < N 3 (\\Du\\ L2{Ql) + \g\ 0]Qs/4 +I[u g>z ,](l) + |m|o;Q 3/4 ) 

+^3 (\fW,Q R (za) + [^71/2,71^3/4) ■ 

To finish the proof of (j2.3|) , it suffices to use Corollary 13.51 

Now we remove the smoothness assumption on the coefficients and data by 
using a standard approximation argument, which we sketched below for the 
completeness. Let A'jjL, m = 1, 2, ... be the mollifications of A a @ . Similarly, 

we define B^ m) , B^ m) , C^ m) , g {m) and / (m) . We know that A^,m = 1, 2, ... 
are uniformly elliptic with the same ellipticity constant v, and as m — > 00, 

{f{m),9(m)) -> (/iff) in L 2 (Qi), 

Moreover, 
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Let V m be the parabolic operator with coefficients A^fL , B? m ) an< ^ C(m) 

instead of A a ^,B a and C. Let v m 6 %|(Qi) be the weak solution to the 
equation 

VmVm = div( 5 - g {m) ) + / - f (m) + ^((Agj - ^ a/3 )^n 

+ (B? m) - B a )u) + (B? m) - B a )D a u + (C (m) - C)u (4.16) 

in Qi with the zero initial-boundary condition v m = on the parabolic 
boundary of Q\. Note that the right-hand side of ()4.16j) goes to zero in L2 as 
m — > 00. It follows from the L2 estimate that ll^mllz^QO + ll-^mllz^Qi) ~~ * ^ 
as m — > 00, which further implies that there is subsequence, which is still 
denoted by {v m }, such that |-Du m | + |v m | — > a.e. in Q\. Let u m = u — v m . 
It is clear that u m is a weak solution to the equation 

V m u m = divc/ (m) + / (m) 

in Q\. By the classical parabolic theory, u m is infinitely differentiable in Q\. 
Thus, by the proof above and by slightly shrinking the domain, we obtain a 
uniform estimate 

\Du m \ . Ql/2 < N(I[u} giZ ,](l) + \g(m)\0;Q 5/6 + l/(m)|o ; Q 6/6 + \\u m IIl 2 (Q 5/6 )) 
< N{I[u g j){l) + |5|o;Qi + |/|0;Qi + \\u m \\ L 2 (Qx)) ■ 

Since Du m — > Du a.e. and u m — > u in L2 in Q±, by taking the limit in the 
above inequality we get (|2.3p . 

Step 2: Continuity of D x /u and U. Next, we prove the second claim of 
the theorem. Fix a point zq G Q1/4 and take 0<r<i?<l/4. 

We define u(R) = u gtZ i{R) + ub,z'{R) + ^A,z'(R) + R J , which is a Dini 
function. By the estimates on u and Du which have already been established, 
we get from (|4.10p that 

/ \D X 'U- (D x iu) Qr[zQ) \ 2 + \U - (U) Qr{zo) \ 2 dz 

< N^r/R)^ f \D x ,u - (D x ,u) Qr{zo) \ 2 + \U - (U) Qr(zo) \ 2 dz 

+ N 4 oj(R)(R/r) d+2 , (4.17) 

where N4 is independent of r, R and zq . Again we take r = tR with r = 
r{d,n,v) sufficiently small such that N\t 21 < 1/2. By an iteration, we 
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obtain from (I4.17D 



/ \ D x'U - (D x ,u) Qr{zo) \ 2 + \U - (U) Qr ( Z0 )\ 2 dz 

<2- h { \D x ,u-(D x ,u) QR{zo) \ 2 + \U-(U) QR{zo) \ 2 dz 



k 

+ N^2-iuj{T k -iR). 
i=i 

We set R = 1/4 and use Lemma I3TT1 to conclude that for r E (0, 1/4) 
\D X 'U - (D x/U ) Qr(zo) \ 2 + \U- (U) Qr{zi)) \ 2 dz < w(r), 



where u> is a Dini function. It then follows from (6.3) of |18j that D x m 
and U are uniformly continuous in Qiu. By using a dilation and covering 
argument, the continuity of D x m and U in Qx/2 follows. 



5. Systems with partially Holder coefficients 

We prove Theorem 12.21 in this section. Thanks to Theorem 12. 1\ we only 
need to estimate the last three terms on the left-hand side of (|2.4p . We 
estimate them separately by using different methods. 

5.1. Estimate of D x /u and U. We fix a point zq £ Q\u and take < r < 
R < 1/4. Following the proof of Theorem 12. 11 we take 71 = 5 and 7 G (5, 1). 
Note that under the conditions of Theorem 12. 2\ we have 

UgA R ) ^ l9]z',8/2,S RS > 

and similar inequalities for ua,z'(R) an d ujb,z'(R)- Owing to (|4.10p . we get 
I \D x ,u - (D x ,u) Qr{zo) \ 2 + \U- (U) Qr{zo) \ 2 dz 
< N 1 (r/R) d+2+ *y [ \D x ,u - (D x ,u) Qr(zo) \ 2 + \U- (U) Qr(zo) \ 2 dz 

+ N 2 R d+2+2S ([g} z , mS + \f\0;Q R (z ) + H5/2,8;Q 3/i + \Du\^Qn{zo)f , 

(5.1) 

8 /2 S 

where N\ only depends on d, n and u, and N2 also depends on the Cj ' 
semi- norms of A and B. 
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Since (15. ip holds for any 0<r<i?<l/4 and 5 < 7, by a well-known 
iteration argument (see e.g., [U Lemma 2.1, p. 86]), 



\D x ,u - (D x ,u) Qr{zo) \ 2 + \U - (U) Qr{z{)) \ 2 dz 

(20) 



/ 

< Nir d+2+26 I lDxiU _ (D x ,u) Ql/4{zo) \ 2 + \U- (U) Ql/i[zo) \ 2 dz 
JQi/a{ z q) 

+ N 2 r d+2+2S ([g] zl ,5/2,6-^ + Ifloi0«(«o) + \u\s/2,6;Q 3/4 + \Du\ . Qr{zo) 



(5.2) 

By the definition of f and U e , we get from (|5.2p . Theorem 12. 1\ Corollary 13.51 
that 

\D x ,u - (D x iu) Qr{zo) \ 2 + \U - (U) Qr{zo) \ 2 dz 

Oo) 

< N 2 r d+2+2S [\g\ z , m5 , Ql + |/|o ;Ql + IM| L2(Ql) ) 2 • (5.3) 

Since (|5.3p holds for any r € (0, 1/4) and zo € Q1/4, by Campanato's char- 
acterization of Holder continuous functions, we obtain 

[D x 'u]s/2,8;Q 1/4 + [U]s/2,5;Q 1/4 < ^2 [\g\z> ,8/2,S;Qi + l/|o ; Qi + IMU a «i)) ■ 
This together with a dilation and covering argument gives 

[D x 'u]s/2,8;Q 1/2 + [U]6/2,8;Q 1/a < N 2 (\g\z' ,8/2,8-^ + l/|o ; Qi + IMU 2 (Qi)) • 

(5.4) 

5.2. Estimate of (u) 1+s . We estimate {u) l+s by modifying the argument 
in [5], which in turn is based on an idea by M. Safonov. The argument in 
[5| uses the maximum principle, which is unavailable for systems. Here we 
use some estimates established in Section [3] instead. 

We denote by Pi the set of all functions p on M. d+1 of the form 

d-l 

p{z) = p(t, x\x d ) = + A 

i=i 

Then we define the first-order partial Taylor's polynomial with respect to 
z' = (t, x') of a function v on at a point z' Q = (to, x' ) as 

d-l 

f} 4 v(z\x d ) := v(z' ,x d ) + ^ A^ ,^)(* 4 ~ 

i=i 

Let 

az') = at,x 1 ,...,x d - 1 )=ri(t+l)l[r 1 (x l ), 

i=l 
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where r\ is a smooth even function on R with a compact support in (—1, 1) 
satisfying 

/ r](t)dt = 1, / t 2 i](t)dt = 0. 

For e > let C e (i, x') = e~ d ~ 1 (^(e~ 2 t, e~ l x') and define a partial mollification 
of v with respect to z' as 

v £ (t, x', x d ) = f f _ i «(s, ?/', x d )C e (t - s, x' - 2/') dy' ds 

v(t — e 2 s, x — ey , x d )((s, y') dy' ds. 

-i 

By virtue of Taylor's formula, it is not hard to prove the following lemma 
for partial mollifications (see, e.g., [HI Chapter 3]). 

Lemma 5.1. Let T G (-00,00] and 5 G [0,2]. Suppose v G cf, /2 ' 5 (R^ +1 ). 
Then for any e > 0, 

e 2 ' 5 sup \D t v £ \ +e 1 ' 5 sup \D x ,v E \ + e 2 ~ s sup \D 2 x ,v £ \ < N[v] z , iS/2 ,s, 

jrpd+l -rrpd+1 ro>d+l 

IKL^-, IKrp M./p 

sup \v - v e \ < N(d,5,'i])e 5 [v] z/) s/2,s, 

where N = N(d, 5, rf) . 

Now we are ready to estimate (u) 1+s . First assume that (|2.1|) holds in 
Rq +1 := {{t,x) G R d+1 : t < 0} and u G cJ +5)/2 ' 1+5 (R^ +1 ). Define 

IMIunif,L 2 , loc := sup ||«||i 3 (Qi(«i))- 

Qi(^i)CK^ +1 

We claim 

Mz',(l+<5)/2,l+<5 < N 2 (\g\ z ' t S/2,5 + \f\o + IM|unif,L 2 , loc ), (5.5) 

where N 2 = N 2 (d,n,u,S, K, [A\ z , jS / 2i s, [B] z ' ,5/2,5) > °- Taker > 0. Let k > 4 
be a number to be chosen later. Denote 

V u = -u t + D a (A a ?(0, x d )D p u). 

Then we have 

V u = div(g + m) + f , 

where 

mi{z) = (A°f(0,x d )-A°f(z))Dp U i, 
3 = (9a), Q a -=g a -B a u, f:=f-B a D a u-Cu. 

Note that 

\m\<N\z'\ s [A] z ,^ 5 \Du\ . (5.6) 
Moreover, since the coefficients of Vo are independent of z' , we have 
V Q u Kr = div(g Kr + rh Kr ) + f r . 
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Let v G T-L^iQnr) be a weak solution to the equation 

f V v = div( - g Kr + m - m Kr ) + f - f r in Q Kr ; 
\ v = on d p Q Kr . 

By Lemma 13.21 we get 

\\v\\ L2iQKr) < iVKr|| -r r +m-m-||L 2 (Q Kr )+ Ar (^) 2 |lf-rilL 2 (Q Kr )- ( 5 - 7 ) 
Corollary 13.51 and a scaling argument yield 

\v\o-,Q Kr/2 < Nkt\q - g Kr + m - m KT \ 0]QKr 

+iV( K r) 2 |f-r|o;Q KI .+iV(Kr)-( d+2 )/ 2 || ? ;|| L2(QKr) . 
This together with (I5.7p . (|5.6p and Lemma 15 . 1 1 yields 

\v\o-,Q Kr/2 < N( K r) 1+s ([Q] z , iS/2)S + [A] z , ms \Du\ + (wrj^lflo). (5.8) 

Let w = u — u Kr — v. Then w satisfies Vqw = in Q Kr . It follows from 
Lemma 13.61 Lemma 13.41 and a scaling that 

\ w ~ ^z' w \o;Qr < ^ 2 H z ',i,2;Q r < N^ 2 \w\ .Q Kr/2 
With the triangle inequality, (15. 8h and Lemma l5.1| we further deduce 

\ W ~ T^w\ .Q r < NK- 2 (\v\ .Q Kr/2 + \U~ U Kr \0;Q Kr/2 ) 

< NK 5 - l r x+s {[Q] z , m5 + [A} z/>s/2>s \Du\o + (Kff-'lflo) 
+ NK 5 - 1 r 1+5 [ U } z , )il+S)/2A+s . (5.9) 

By Lemma [53] we also get 

- f^ r \ , Qr < Nr\[Dl,u^] . Qr + [A^1o ; qJ (5.10) 

< A^-V +, V]^ (1+a)A1+a . 

Take p = f l z ,w + f l z ,u Kr G Pi. Then from ipTB]> . (JESJ) and (IBTTOl) . we get 

|« - P|o ; Q r < Mo ; Q r + K r - T*/ U Kr |o ; Q r + |w - r^tw| 0; Q P 

< iV(Kr) 1+5 ([ 5 ] zV/2i(5 + [A] z , ms \Du\ Q + {nr) 1 - 5 ^) 
+ N K 5 - 1 r 1+5 [u] z , Xl+5) ^ 1+s . 

By a shift of the coordinates, we have for any zq G Kq +1 and r > 0, 
inf \u - p\ . Qr{zo) 

pGPl 

<^ 1+5 ([g] yi5A5 +[A],, iV2i5 |^|o+(Kr) 1 - 5 |f|o)+A r iK 5 - 1 [«] z ',(i+5)/2,i+5, 

(5.11) 
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where Ni = Ni(d,n,p) > 0. We fix a large k = K{d,n,v) such that 
Nik 5 ^ 1 < 1/2. For any r < 1/k, we get from (15. lip 

r^ s inf \u - p\ 0;Qr ( Z0 ) 

pGPl 

< ^([q]z',5/2,5 + [A] z > t S/2,s\ D u\o + Iflo) + l/2[«]y,(l+<S)/2,l+<S- 
By the definition of f and 0, 

r- 1 - 6 inf \u - p\ . Qr{zo) 
pePi 

< N ([g]z>,6/2,5 + Mo + Hz',5/2,6 + \ Du \o + |/|o) + 1/2 ,(l+,5)/2,l+ ( 5 

< AfcflffU'.tf^a + |/|o + ||«||unif,i 2) i oc ) + l/ 2 M z / i (i +5 )/2,l+5 ) ( 5 - 12 ) 

where = N2(d,n,i',8,K, [A] z / t s/2,Si [B] z ', 8/2,5) > 0- I n the last inequality, 
we used Corollary 13.51 and Theorem 12.11 On the other hand, for r > 1/k, 
clearly we have 

r_1 ~ 5 ' m ! \ u -p\o;Q r (z ) < N\u\q. 
pGPl 

Thus, ([5.12p holds in any case. 

Now by first taking the supremum in (|5.12p with respect to r and zq, and 
then using the equivalence of parabolic Holder semi-norms similar to []Tj 
Theorem 8.5.2], we obtain 

Mz',(l+<5)/2,l+<5 < N 2(\g\z>, 5/2,5 + l/lo + |M|unif,L 2 ,ioc ) + l/ 2 [ n ]z',(l+5)/2,l+5' 

which implies (|5.5p under the assumption that u € C^, 1 1+<5 (Rq +1 ). 
Next, we localize the estimate. Let u be a weak solution to (|2.ip in Q\ 

and assume u € C^ +S ^ 2 ' 1+S (Q3/4). We take a cutoff function <j such that 

q = in Mq +1 \ Q3/4 and s = 1 in Qi/2- Then u = u(, satisfies 

Pu = div<? + / inK^ +1 , 

where 

f = fq - - + (A Q/3 ^n + (B a + B a )u) D a s. 

We then have (|5.5p with u, f, g in place of u, f, g. This together with Theo- 
rem 12.11 and Corollary 13.51 yields 

[ u )z',(l+S)/2,l+S;Q 1/2 < N 2 [\g\z',S/2,S;Qi + l/|o;Qi + \\ u \\l 2 (Qi)) ■ ( 5 - 13 ) 

Finally, we drop the assumption that u G C^} +S ^ 2 ' 1+S {Q3/4) by using the 
same approximation argument as in the proof of Theorem 12.11 Combining 
(15. 4p and (15.13p . the theorem is proved. 
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6. Non-divergence form equations 

In this section, we are concerned with non-divergence scalar parabolic 
equations 

Pu := -u t + a aP D af) u + b a D a u + cu = f, (6.1) 

The coefficients a a/3 , b a , and c are assumed to be measurable and bounded 
by K, and the leading coefficients a a @ are uniformly elliptic with ellipticity 
constant v: 

HCI 2 < a afS Cf, \a al3 \ < v~ l 
for any £ G M. d . The aim here is to prove estimates in the same spirit as 
in Theorems 12.11 and 12.21 for non-divergence form equations. We state the 
main results of this section as follows. Theorem 16.21 improves Theorem 2.14 
[5] and Theorem 3.1 [19J in the case q = d — 1. 

Theorem 6.1. Let a,b,c G C$ ini and f G C$ ini (Qi). Assume that u G 
W2' 2 (Qi) is a strong solution to f)6. 1 1) in Q\. Then we have 

Ml,2;Q 1/2 < N(I[Uf^](l) + |/| 0;Ql + \\U\\ L2{Q1) ), (6.2) 

where N = N(d,v,K,u atZ r,ujb :Z >,Lj C:Z >). Moreover, ut and D xx /u are contin- 
uous in Q\/2- 

Theorem 6.2. Let 5 G (0,1), a, 6, c G C S J 2 ' 5 and f G cf, /V (Qi). Assume 
that u G W 2 ' (Qi) is a strong solution to (|6.ip in Qi. T/ien we /lave 

|«|l,2;Q 1/2 + ['W i ]5/2, ( 5;Q 1/2 + [D xx 'u] S /2,5;Q 1/2 < N (\f\ z < j /2,5;Qi + IMIl 2 (<2i))> 

where N = N(d,5, u,K, [o] z / j5 / 2) a, [6] 2 ',5/2,<5, [c] «/, ,5/2,5 )■ 

Remark 6.3. From Theorems 16. II and !6.2l we also obtain the corresponding 
results for non-divergence form elliptic equations, as in Corollaries 12.61 and 
12^1 

Remark 6.4. We consider the case when the coefficients and data are piece- 
wise Holder continuous. Suppose that Q\ is divided into M laminate sub- 
domains T>\ , T>2, ■ ■ ■ , Dm by M — 1 parallel hyperplanes with the common 
normal direction (0, • • • ,0,1). Assume a,b,c and / are in C s ' 2 ' S (T>i) for 
each i = 1, 2, • • • , M, but may have jump discontinuities across these hyper- 
planes. By Theorem 16.21 and a covering argument, for any e G (0,1/2), ut 
and D TT iu are in C s/2fi {Q\- e ). Moreover, by the same theorem, u and Du 
are m C s/2 ' 5 (Qi- £ ). Now restricted to each V { n Qi- £ , since 

D dd u = {a dd )- 1 If + iH- b a D a u -cu- a aP D aP u 

we conclude that D dd u G C s/2 ' 5 (T>i n Qi- £ ), and hence u G C 1+5/2 ' 2+5 (Vi n 
Qis). It is worth noting that the C 1+s / 2 > 2+s norm of in each sub-domain 
T>i n Qis is independent the number M. 
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The proofs of Theorems 16.11 and 16.21 follow the line of the proofs of The- 
orems 12.11 and 12.21 an d in fact are simpler. First, by using the argument 
in Section we obtain from the main result of [TO] the following interior 
estimates, which is analogous to Corollary! 



Lemma 6.5. Let q € [2, oo). Assume a al3 are VMO in z' , u € C^ c 
satisfies (16. ip in Q\, where f € L q (Qi). Then there exists a constant 
N = N(d, v, K, oj a ,z' > q) such that 

l|u|l ^' 2 (Qi/ 2 ) - N (W u Wl 2 (Qi) + ll/llr,(Qi))- 
In particular, if q > d + 2, it holds that 

M(l+ 7 )/2,l+ 7 ;Q 1/2 < N{\\u\\ L2{Ql) + 

where j = 1 — (d + 2)/q. 
We denote 

P (u) = -ut + a aP D a pu. (6.3) 
The next lemma is an analogy of Lemma 13.61 

Lemma 6.6. Let 7 G (0,1). Assume u E C^ c satisfies Pqu = fo in Q\, 

where f$ and a a @ are independent of t and x' . Then there exists a constant 
N = N(d,j,v) such that 

M(i+7)/2,i+7;Q 1/2 ^IMImqi)' ( 6 - 4 ) 

[AH(l+7)/2,l+7;Qi/2 ^ N \\ D xx'U\\L 2 {Q x y (6.5) 
Proof. Note that in Q\ 

Po{DiD j xl u) = 0, 

wherever i + j > 1. Then (16. 4p follows from Lemma 16.51 applied to Uf. For 
any 1/2 < r < R < 1, applying the same lemma to D x m — (D x /u)q r gives 

[A^«](i+7)/2,l+7;Qr < N\\D x ,u - {D x ,u) Qr \\ L2{Qr) . (6.6) 

To bound the right-hand side of (|6.6p . we use an idea in [6] to utilize the 
divergence structure of the equation after making a change of variables. Let 

fX d I 

y d = ip{x dy = —-ds, yP=xP, /3 = 1,--- ,d-l. 
J a ad (s) 

It is easy to see that ip is a bi-Lipschitz function and 

5 < y d /x d < r 1 , D y d = a dd (x d )D xd . 

Denote 

v(t,y',y d ) = u(t,y',<p- l (y d )), a«(y d ) = a^VV)), 

/V) =/(<p- 1 (v d ))- 

Define a divergence form operator Pq by 

P 0V = - Vt + D d ( w D d v)+J2D P i^ d —D d v)+ D a (d^D,v). 
V ° J p=l V a J a ,p=l 
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Clearly, Pq is uniformly nondegenerate and v satisfies Pqv = f in some 
stretched cylindrical domain. Since Po(D y /v) = 0, one can use Lemma 13.31 
i) to estimate the mean oscillation of D y /v in each parabolic cylinder by the 
integral of D yy /v in the same cylinder. To finish the proof of (|6.5|) . it suffices 
to use a covering argument bearing in mind that the metrics in ^-coordinate 
and y-coordinate are comparable. □ 



Now we are ready to prove Theorems 16.11 and I6.2L 

Proof of theorem 1 6. 1\ First we assume that b = c = 0. Also, arguing as 
before we may assume that o are infinitely differentiable and u has bounded 
derivative up to fourth order in Q 3 / 4 . We take < 7 < 1. Fix a point 
z o € Q3/4, and take < r < R < (3/4 - |x |)/4. Now take z[ € Q' R {z' ) and 
denote 

P z ,u = -u t + A a P(z' 1 ,x d )D a pu. 

Then we have 

p z[ u = f + m, 

where 

m{z) = {a af3 {z' llX d ) - a al3 (z))D a pu. 
Let v be the strong solution to the equation 

f P z[ v = f(z) - f(z[,x d ) + m(z) in Q 2 r(zo); 

[ v = on d p Q 2 R{z ). 

We rewrite P z j as a divergence form operator as in Lemma 16.61 then use 
Lemma 3.2 to get 

1 1 V 1 1 L 2 (Q 2R (z ))< NR2 1 1 / " / (*i ,X d ) + m 1 1 L 2 (Q 2R (zo))- 

This together with Lemma 16.51 gives 

II^HmQkOso)) + W v t\\L 2 (Q R (z )) <N\\f- f(z[,x d ) + m\\ L2iQ2R{zo)) . (6.7) 

Let w = u — v. Then w satisfies P z ^w = f(z' 1 ,x d ) in Qr(zo). It follows from 
Lemma 16.61 and a suitable scaling that 



/ 

JO,., 



\D xx ,w - {D xxl w) Qr{za) \ 2 + \ w t - (w t ) Qr(zo) \ 2 dz 

Oo) 



<N{r/R) d+2+2 ^ f \D xx ,w\ 2 + \w t \ 2 dz. (6.8) 
•>Qr(zq) 

Define 

Hx d ) = £ j\a dd {z[,r))- 1 (a*? \z[,t)(D^w) Qr(zo) - (wt) Qa(zo) ) drds, 



and 



w := w - t(wt) QR{zo) - -x a x^{D olP w) QR{zo) + h(x d ). 
Then, we have 

D xxl w = D xx >w - (D xx >w) Qr ( Zo ), w t = w t - {wt)Q R ( Z0 )- 
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It is easily seen that w also satisfies P Z 'W = in Qr(zq). We substitute w 
in (|6.8p by w to get 



/ \D xx ,w - (D xx ,w) Qr ( Zo) \ 2 + \w t - (w t ) Qr ( Z0 )\ 2 dz 

J Qr(zo) 



< N{r/R) d+2+2 ^ / \D xx ,w - (D xx ,w) Qr(zo) \ 2 + \w t - K) Qfl(zo) | 2 tte. 

(6.9) 

We combine (|6.7p with (|6.9p and use the triangle inequality to obtain 

/ \Dxx'U- (D xx ,u) Qr{zo) \ 2 + \D t u - (D t u) Qr{zo) \ 2 dz 

J Qr(zo) 

< N l {r/R) d+2+2 ^ I \D xx ,u - (D xx ,u) Qr{zo) \ 2 + \D t u - (D t u) Qr{zo) \ 2 dz 

J Qr(zo) 

+ N\\f-f(z' 1 ,x d )+m\\l 2{Q2R{zo)) , (6.10) 

where N\ = N\(d,v). Now we take average of both sides of (|6.10p with 
respect to z[ € Q'r(z' ) to get 



/ \D xx >u - {D xx m) Qr ( ZQ) \ 2 + \D t u - (D t u) Qr(zo) \ 2 dz 

JQr(zo) 

< N x {r/R) d ^ + ^ [ \D xx ,u - {D xx ,u) Qr{zo) \ 2 + \D t u - {D t u) Qr{zo) \ 2 dz 
+ NR d+2 ujl z ,(2R) + NR d+2 u J l z ,(2R)\D 2 u\l Q2R{zo) , (6.11) 



As in Section we immediately get (|6.2p from (|6.1ip by using an iteration 
argument and Lemma 16.51 Then the argument in Step 2 of the proof of 
Theorem 12.11 shows the continuity of D xx /u and Ut- 

In the general case, we move all the lower order terms to the right-hand 
side: 

— ut + a a/3 D a /3U = f — b a D a u — cu := f. 
From the proof above, we get 

Ml, 2 ;Q 1/2 < N{I[U^]{1) + IflojQi + \\u\\ L , 2{Ql )). 

To bound the first two terms on the right-hand side, it suffices to use Lemma 
16.51 and the assumptions on b and c. □ 



Proof of Theorem \6.2l First we assume that b = c = 0. We fix a point 
z ^ Qi/4 an d take < r < R < 1/8. Following the proof of Theorem 16. 1( 
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we take 7 £ (5, 1). Owing to (|6.10p . we get 

f in ... /-n „A |2 , |„. /„, ^ |2 



[Am/U - (I> a!a! '«)Q r ( jeo ) | 2 + \u t ~ (ut)Q r (z ) ? dz 



< N 1 (r/R) d+2+ ^ / \D xx ,u - (D xx ,u) Qr[zo) \ 2 + \u t - (u t ) QR{zo) \ 2 dz 

JQr{zq) 

+ N 2 ([fU, S /2,5;Qi + \D 2 uW,Q 2R { Z0 )) 2 R d+2+2& 1 (6-12) 

S /2 S 

where N\ only depends on d and is, and .A/2 also depends on the Cj ' semi- 
norm of a. Since (|6.12p holds for any 0<r<i?<l/8 and 5 < 7, by a 
well-known iteration argument (see e.g., [9l Lemma 2.1, p. 86]), 

IAkb'W - (A^'W^Oo) | 2 + \u t - (u t ) Qr ( Z0 ) | 2 dz 

Qr(zo) 

< N ir d+2+2S [ \D xx ,u - (D xx ,u) Qi/dzo) \ 2 + \u t - (u t ) Ql/4(zo) | 2 dz 

+ N 2 r d+2+25 {{fU, s/2j5 . Ql + |^ 2 n| ;Q 2B ( Z0 )) 2 • (6.13) 
We get from (|6.13j) . Theorem 16.11 and Lemma 16.51 that 



/ \D xx ,u- (D xx ,u) Qr ( Zo) \ 2 + \ut - {ut)Q r ( ZQ )\ 2 dz 

JQr(z ) 

< N 2 r d+2+2S (\f\ z , mS;Ql + \\u\\ L2(Ql) ) 2 . (6.14) 

Since (|6.14p holds for any r € (0,1/8) and zq € Q1/4, by Campanato's 
characterization of Holder continuous functions, we obtain 

[D xx >u] S / 2t S;Q 1/4 + [ut\5/2,S;Qi/4 - N<2 {\f\z',S/2,6;Qi + IWIUa(Qi)) ■ 
This together with a dilation and covering argument gives 

[D xx 'U] S / 2 ,S;Q 1/2 + [ u t]s/2,S;Q 1/2 < ^2 (\f\z',S/2,S;Qi + IMU 2 (Ql)) ' 

In the general case, we move all the lower order terms to the right-hand 
side and ague as before. □ 

We finish this section by proving the following "partial" Schauder esti- 
mates. These results generalize Theorem 5.1 [5 J and Theorem 2.1 [19J for 
the Poisson equation. 

Let us introduce a few more notation. Let q be an integer such that 1 < 
q < d— 1. We distinguish the first q coordinates of x from the rest and write 
x = (x, x), where x = (x l , ■ ■ ■ , x q ) and x = (x q+1 , • • • , x d ). We also denote 
z = (t,x). As in Section [21 we introduce the partial Dini continuous space 
CP ini (£>) and the partial Holder space C~^ 2 ' S , as well as their corresponding 
norms. 
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Theorem 6.7. Let a = a(x d ) be a measurable function of x d alone. Let 
b,c G CP ini and f G CP im (Qi). Assume that u G W 2 1,2 (Qi) is a strong 
solution to (|6.1|) in Q\. Then we have 

M0;Q 1/2 + \D xx u\ -Q 1/2 < N{L[u ftZ i}{\) + 1/lojQi + ||u|| L2 ( Ql) ), 

where N = N(d,iy,K,u}b tZ ,u} CtZ ). Moreover, ut and D xx u are continuous in 

Ql/2- 

Theorem 6.8. Let 5 G (0,1), a = a{x d ) be a measurable function of x d 

alone. Let b,c G C S / 2 ' S and f G C S / 2 ' S (Qi). Assume that u G W 2 1,2 (Qi) is a 
strong solution to (|6.ip in Q\. Then we have 

Wt\8/2,S;Q 1/2 + \D xx u\s/2,8;Q 1/2 < ^ (\f\^ 5 /2,8;Q 1 + IMIl 2 (Qi))> 

where N = N(d,5, u,K, [b] iiS / 2 ,s, [c]z,5/2,s)- 

For the proofs of Theorems 16.71 and 16.81 first we note that similar to 
Lemma 16.61 if /o is independent of z, then (|6.4p still holds and we have 

[Ac u ](i+7)/2,i+7;Qi/2 < N\\D xx u\\ L2 ( Ql) . 

Following the proof of Theorem 16.11 first we assume there is no lower order 
terms. Let v be the solution of 

j P v = f(z) - f{zi,x) mQ 2 R(z ); 

| v = on d p Q 2 R(z ), 

where Pq is defined in (|6.3p . Then as before we have 

W d2v \\l 2 (Q r (z )) + IMUaWafr,)) < ^11/ - f(zi^)h 2 (Q 2R (z ))- (6-15) 
Moreover, w = u — v satisfies 

\D xx w - {D xx w) Qr{zo) \ 2 + \w t - (w t ) Qr ( Z0 )\ 2 dz 

r(z ) 

< N(r/R) d+2+ ^ [ \D xx w - (D xx w) Qr(zo) \ 2 + \w t - (w t ) QR{zo) \ 2 dz. 

(6.16) 

Combining f)6. 15|) with (|6,16p and taking average in z\ gives 
\D xx u - {D xx u) Qr[za) \ 2 + \D t u - {D t u) Qr{zo) \ 2 dz 

-(zo) 

< N x {r/R) A+2 ^ f \D xx u - (D xx u) Qr(zo) \ 2 + \D t u - {D t u) Qr{zo) \ 2 dz 
JQr{zo) 

+ Noo 2 f:Z (2R)R d+2 . 

Now a standard iteration argument finishes the proof. 

We remark in Theorems 16.71 and 16 . 81 if b, c and / are assumed to be regular 
only with respect to x, we can still get the estimate of D xx u by dropping 
the ut terms in (|6.15p and (|6.16p and replacing f(zi,x) by f(t,X\,x). 
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7. Appendix 

In the appendix, we give a proof of Lemma 13.41 Let Ao be the constant 
in Theorem 2.2 of [4j. Let 

r k = l-2~ k , Q k = (-r 2 k ,Q) x B rk , k = l,2,---. 

Then we find ( k {t,x) G C^(R d+1 ) such that 

f 1 on Q k 

[0 on R d+1 \(-rl +l ,rl +1 )xB rk+1 



and 



\D( k \ < N2 k , \(C k ) t \ < N2 2k , \D\ k \ < N2 2k . 
Observe that, for X k > Xo, 

(V - X k )(( k u) = div g k + f k in R$ +1 , 

where 

d 



9k = (SfcJ , 9k a = Ck9a + ^ A a ^uD^( k , 



13=1 

d d 



fk = Ckf 9aD a (k + Aaf>D p uD aC,k - uD t ( k - X k uC, k . 

a=l a,P=l 

Then by Theorem 2.2 of 

W D (Cft«) Hl,(R*W) ^ N (ll^ll Lp (M d + 1 ) + ^^IIMIl, (<+')) 

< iV (2 k + A^V* + Af ) ||u|| M<w + iVA^II/lli,^, 



+ AT (l + A, 1/2 2 fe ) || 5 || MQl) + NX k 1/2 2 k \\D (Oh-iu) H^i). 

Set 

A = H-D (Cfc«) ll Lp(R d+i } , = IMIi^qo, = II^IU P (QO' T = ll/lli»(Qi)- 
Then 

A<ivf2 fe + A- 1/2 2 2fe + AfU 



+ N ( 1 + A, 1/2 2 fc ) + iVA" 1/2 F + .ZVA7 1/2 2 fe A +1 . 



By multiplying e both sides and summing up with respect to k, we have 

oo oo 

Y^e k A k <NBJ2 (2" + A^ 1/2 2 2fc + Af ) e k 
k=i k=t 

oo oo oo 

+ NQ <T (l + A fc - 1/2 2 fc ) e fc + A fc " V2 ^ + M E A- 1/2 (2 £ )*A+i- 

fc=l fc=i fc=i 
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1/2 

We may certainly assume N\ > A . Now set 

e = 1/4, A*/ 2 = N 1 2 k+2 . 

Then 

oo oo 

jr (2 fe + A, 1/2 2 2fc + Af ) e fc < N, £ (l + A^ 1/2 2 fc ) e fc < AT, 
fe=i fc=i 

oo oo oo oo 

^A fc 1/2 e fc < iV, iV 1 ^A- 1 / 2 (2e) fe A + i = E £fe+1 ^+i = 

fc=l fc=l fc=l fc=2 

Therefore, 

oo oo 

J> fc A <N(B + G + F) + ^2e k A k . (7.1) 

fc=l fc=2 

On the other hand the summations above are finite because 

A k <N2 k \\u\\ Lp(Ql) +N\\Du\\ Lp(Ql) . 
Then the inequality (|7.ip implies that 

\\Du\\ Lp{Qi/2) <Ar< N(B+g+J r ) = N(\\u\\ LAQl) + \\g\\ LpiQl) + \\f\\ Lp{Ql) ). 

Finally, the estimate of ||' u t||ji- 1 (Q 1 2 ) follows from the above estimate and 
the equation f)2. 1 1) itself. The lemma is proved. 
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